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A way to efficiently compute helicity amplitudes for arbitrary tree-level scattering processes in QCD is presented. 
The scattering amplitude is evaluated recursively through a set of Dyson-Schwinger equations. The computational 
cost of this algorithm grows asymptotically as 3", where n is the number of particles involved in the process, 
compared to n! in the traditional Feynman graphs approach. Unitary gauge is used and mass effects are available 
as well. Additionally, the color and helicity structures are appropriately transformed so the usual summation is 
replaced by the Monte Carlo techniques. 



1. INTRODUCTION 

QCD processes with many external legs are of 
much interest, both for testing QCD in different 
settings and as backgrounds for new physics pro- 
cesses at the Fermilab TeVatron and at the CERN 
LHC. However, the estimation of multi-jet pro- 
duction cross sections as well as their character- 
istic distributions is a difficult task. Perturbation 
theory based on Feynman graphs runs into com- 
putational problems, since the number of graphs 
contributing to the amplitude grows like n\. The 
counting of the graphs itself becomes a problem 
let alone their evaluation and the computation 
of the color and helicity structures which is an 
additional source of computational inefficiencies. 
Over the last few years, new recursive algorithms 
based on Dyson-Schwinger equations |1I2I3| or 
on field equations 6 7| have been developed 
in order to overcome the computational obsta- 
cles. Very recently also on shell recursive equa- 
tions have been proposed |8I9| . The amplitude 
calculated using Dyson-Schwinger recursive equa- 
tions, which avoid Feynman diagrams, results in 
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a computational cost growing asymptotically as 
3". Here, the off-shell subamplitudes are intro- 
duced which are combinations of parts of Feyn- 
man graphs. For those subamplitudes a recursion 
relation has been obtained which enables to ex- 
press an n-particle amplitude in terms of all sub- 
amplitudes, starting from 1—, 2—, ... up to (ji—l) 
particles. The color and helicity structures cor- 
responding to each subamplitude have a simpler 
form. Moreover, they can be appropriately trans- 
formed so the usual summation can be replaced 
by the Monte Carlo one. 

In this article the algorithm based on Dyson- 
Schwinger recursive equations is presented and 
used in order to efficiently obtain cross sections 
for arbitrary multi-jet processes. 



2. DESCRIPTION 
RITHM 



OF THE ALGO- 



To illustrate how the algorithm works let 
us present a simple example, the gluon self- 
interaction: 



(1) 



Here pi,p2,P3,Pi represent the external momenta 
involved in the scattering process, taken to be in- 
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coining. The subamplitude with an off-shell gluon 
of momentum P has contributions from three- 
and four-gluon vertices only. To reduce the com- 
putational complexity down to an asymptotic 3", 
we replaced the four-gluon vertex with a three- 
gluon vertex by introducing an auxiliary field rep- 
resented by the antisymmetric tensor _ff^'^. The 
recursion for the gluons now changes only slightly. 
However, we have an additional equation for the 
auxiliary field. Diagrammatically, the full con- 




Figure 1. Recursion equations for gluon self- 
interactions. The auxiliary field is used to replace 
the four-gluon vertex in order to reduce computa- 
tional complexity of calculations. 



tent of Dyson-Schwinger equations is presented 
in Fig. n and the corresponding equations have 
the following form: 

P=Pl+P2 

+ ^ E KxMAc(Pl)H^'BiP2) 

P=Pl+P2 

-Hac iP2 ) ^CB (Pi ) } e (Pl , P2 ) 

P=Pl+P2 



{A'AciPi)A'cBiP2) - A^c(P2)A^b(pi)} 

where Xf^^^P is the new auxiliary field-gluon- 
gluon vertex: 

and A,B,C — 1,2,3, and e{pi,p2) is the Fermi 
sign factor. In order to be more transparent we 
will write explicitly the three-gluon vertex part 
as well as the auxiliary one, suppressing the color 
indices and using the light-cone representation: 

A^{P) ^ {A{p,) ■ A{p2)){p2 - PiY + {Pi ■ A{P2)+ 
P ■ A{p2))AP{p,) - {p2 ■ + P ■ Aip,))A^{p2) 
+A,{pi)H'^"{p2) - AMH'"'iP2) 
H^%P) A^{pM''{P2) - A^{pi)A^^{j>2) (3) 

where all momenta are taken to be incoming. 
These equations, the off-shell fields, are the 
main building blocks of the gluon self-interaction, 
which is to be constructed iteratively. After n — 1 
steps, where n is the number of particles under 
consideration, one can get the total amplitude 

A{pi,P2, ...^Pn) = AAB{Pi) ■ AAB[Pi), (4) 

where 

P^ = Y.P^ (5) 

so that Pi + Pi ~ 0. The hat denotes functions 
given by the expression from the previous step ex- 
cept for the propagator term. This is because the 
outgoing momentum Pi must be on shell and the 
propagator term is removed by the amputation 
procedure. The iteration begins with the initial 
condition for the external particles. For gluons 
we have 

^ABiPi) = £tiPi)^AcSDB, (6) 

where i enumerates the external particles, i = 
1, n, e denotes the polarization vector and A = 
±1. 

In order to label and systematically control 
the momenta of external particles and their rel- 
evant intermediate combination, we give all mo- 
menta in the binary representation, see e.g. 
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For the process under consideration, and the ex- 



ternal particles with momenta 



1,2,3,4 



we assign to the momentum a binary vector 
m = (mi, 777.2, TO3, 7714) with components which 
are either or 1 as follows: 



4 

P'^^Y.m.pf. 



(7) 



The binary vector can now be uniquely repre- 
sented by the integer 

4 

777 = ^2^-1777,. (8) 

i=l 

In particular, one can write: 
ff(l)5(2)^g(4)5(8). 

All subamplitudes can now be replaced by the 
corresponding integers 

A^(P) ^'^ (777), 

This representation allows us to establish a natu- 
ral ordering of the momenta based on the notion 
of level, defined simply as 



I = y^mi 



(9) 



All external momenta are at the level 1, whereas 
the total amplitude corresponds to the unique 
level 77 = 4. With levels we can see the natural 
iteration path of the equations. We start from 
the subamplitudes at the level 1 which are the 
external momenta together with the initial con- 
ditions, then go to the subamplitudes at the level 
2 and so on until we reach level 77 = 4. The Fermi 
sign factor is also expressed in the integer number 
language 



£(^"1,^2) 6(7771,7712) 

and we use the following formula 

6(7771, 7772) = (-1)X(™1.™2) 

with 

2 i-l 

X(777i, 7772) = X! X! 

i—n J — 1 



(10) 

(11) 

(12) 



where the hat means that this particular compo- 
nent is equal to if the corresponding external 
particle is a boson. This sign factor takes into ac- 
count the sign change when two identical fermions 
are interchanged. 

Contrary to orginal HELAC |1I2| . the computa- 
tional part consists of only one step, where cou- 
plings allowed by the lagrangian defined by fusion 
rules are only explored, see Ref. ^U] for technical 
details. Subsequently, the helicity configurations 
are set up. There are two possibilities, either ex- 
act summation over all 2"^ x 2"^ x 2"^ helicity 
configurations, where 711 is the total number of 
quarks, 772, the total number of antiquarks and 
773, the total number of gluons; or Monte Carlo 
summation. For example for a gluon the second 
option is achieved by introducing the polarization 
vector 



e^(p) = e'^4(p)-He-'V(p), 



(13) 



where (f) G (0, 27r) is a random number. By inte- 
grating over 4> we can obtain the sum over helic- 
ities 



1 



2tt 



^ £^(p)(£^(p))*. 
A=± 



The same idea can be applied to the helicity of 
quarks and antiquarks. 

Finally, the color factor is evaluated iteratively. 
Once again, we have two options. Either we pro- 
ceed by computing all 3"' x 3"« color configu- 
rations, where the gluon is treated as a quark- 
antiquark pair and 77^, 77^ is the number of quarks 
and antiquarks respectively, or particular config- 
urations are chosen by the Monte Carlo method 
(for more details see Ref. prO) 1 . 

For the process g{Pi)g{P2) 9{Pi)g{Ps) 
which we are using through this section as an ex- 
ample the momenta involved in the various levels 
of calculation are the following: 

1st Level: 

Pi = (0001) = pi, P2 = (0010) = P2 
Pi = (0100) = P3, Ps = (1000) = Pi 
2nd Level: 
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P6 = (0110)-P2+P3 
PlO = (1010) =P2+P4 
Pl2 = (1100) =P3+P4 

3th Level: = (1110) = P2 + + P4 

4th Level: P15 = (1111) = pi + P2 + P3 + P4- 

Note that we have chosen the particle number 
1 as our ending point and, therefore, we have 
only computed all the relevant momentum com- 
binations where the momentum pi does not ap- 
pear. This excludes all odd integers between 1 
and 2" — 2. The momentum pi is combined only 
once, in the last n— th level. Starting from the 
second level we get: 

H^'^ie) ^ X^"'P''Ap{2)A^{4) 
yl^(lO) - V''P''ilO,2,8)Ap{2)A^{8) 
H^'^ilO) - X^"'P''Ap{2)A^{8) 
yl^(12) - V''P''il2,4,8)Ap{4)A„{8) 
H^^iU) - XP''P''Ap{4)A„{8) 
At level 3: 

A^iU) - F^'"'(14,2,12)Ap(2)A^(12)- 
XP-P-A,{2)Hp,{12) + VP''^U,4,10)Ap{4)A,{10) 
-X^^''"A,(4)iJp,(10) + yp'^iU, 8, 6)Ap(8)A,(6) 
~X^^•'p-A,,{8)Hp,{6). 

At the last, 4-th level, the total amplitude is com- 
puted by combining this subamplitude with the 
remaining one, which describes the momentum pi 
and is simply given by 

^(15)~ A^(l)-A''(14). 



Table 1 

Results for the total cross section from 4 up to 8 
gluons, CTt corresponds to summation over all pos- 
sible color configurations, while a!^'^ corresponds 
to Monte Carlo summation. 



Process 


(It ± error (nb) 


o-"*^ ± error (nb) 


99 ^ 2g 


4611.55 ± 38.13 


4627.18 ± 33.28 


99 ^ 3g 


152.444 ± 2.490 


152.137 ± 2.822 


99 ^ 4g 


12.9072 ± 0.4070 


12.6137 ± 0.4619 


99 5g 


1.04254 ± 0.05300 


1.04446 ± 0.10390 


99 ^ 6.9 


0.07577 ± 0.00597 


0.07261 ± 0.00516 



3. NUMERICAL RESULTS 

As an example the algorithm has been used to 
compute total cross sections for multiple jets pro- 
duction. The CMS energy was chosen y/s = 14 
TeV and the following cuts were applied 

PT, > 60 GeV, \ri,\ < 2.5, AR > 1.0 (14) 



where px — yPx + Py the transverse momen- 
tum of a jet, 7] — — lntan(0/2) is the pseudo- 
rapidity. Many methods can be used to define 
what is meant by a jet of hadron. One commonly 
used is the 'cone' description of a jet which is the 
transverse energy, Et, concentration in a cone of 
radius 



AP = 
with 



A ^ f P^^P^j + PviPvj \ 
A$,;, = arccos ^— ^ . 

V PTiPT, ) 

All results are obtained with a fixed strong 
coupling constant calculated at the Mz scale. 
There are several parameterizations for the par- 
ton structure functions, we used CTEQ6 PDF's 
parametrization 11 12^. For the phase space gen- 
eration we used either PHEGAS ^Hl or a flat phase- 
space generator RAMBO 14 . 

The results for the total cross section from 4 up 
to 8 gluons are listed above in Tabled We give 
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applied. 

The next result we present is an example of 
summation over all color configurations in com- 
parison to Monte Carlo summation. In Fig. [3 the 
X-axis variable is 

all 

z = |Af,|V^|M,p 

i 

where | Af / 1 ^ is the square matrix element for one 
particular color configuration normalized to the 
sum of all possible is plotted. The agreement is 
easy visible. 

4. SUMMARY 

An efficient tool for automatic computation of 
helicity amplitudes and cross sections for multi- 
particle final states in QCD has been presented. 
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Figure 2. Distribution in z = |M/p/^°" |M,p, 
where |M/p is the square matrix element for one 
particular color configuration normalized to the 
sum of all possible is plotted. The upper plot cor- 
responds to the gg —>■ gg process, the lower one to 
the gg ggg. Solid line crosses denote summa- 
tion over all color configurations whereas dashed, 
the Monte Carlo summation. 



the 'exact' result with summation over all pos- 
sible color configurations (10^ generated events), 
(7t as well as the result obtained with Monte Carlo 
summation over color (4 x 10^ generated events) 
ct"°. In both cases a Monte Carlo over helicity is 
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